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Abstract 

0^ , An ultracold wave packet of bosonic atoms loaded into a travelling laser wave may form a many- 

0^ ■ atom soliton. This is disturbed by a homogeneous force field, for example by the inevitable gravi- 

' tation. The wave packet is accelerated and therefore the laser frequency appears to be chirped in 

^ |. the rest frame of the atoms. We derive the effective nonlinear Schrodinger equation. It shows a time 

D ' dependent nonlinearity coefficient which amounts to a damping or antidamping, respectively. The 

^0 ' accelerated packet solution remains a soliton which changes its shape adiabatically. Similarly, an ac- 

^y-j , tive shaping can be obtained in the force-free case by chirping the laser frequency thus representing 

a way of coherent control of the soliton form. The experimental consequences are discussed. 
y—i \ PACS numbers: 42.50.Vk; 32.80.-t; 03.75.Be. 

^ . to be published in Europhys. Lett. 

oo ; 

' After the recent realizations of Bose-Einstein condensation in ultracold and dilute gases [1-4] , it may 

, be expected for the near future that a new source of atoms will be available which provides atoms of 

' high flux and coherence. The output coupler of such an atom laser is based on atom-light interaction 

, : Coherent coupling of atoms in two distinct hyperfine states by microwaves is used in |^ to transfer 

■ the Bose condensate from a trapped to an untrapped state. A significant consequence of the quantum 

field theory of ultracold bosons interacting with photons j^, ||] is that two-body interactions like the 
dipole-dipole interaction lead to nonlinear effects in the propagation of wave packets. This opens a new 
exciting field of research in atom optics called nonlinear atom optics. It has been worked out by several 
'"^ ' authors |ri| that due to this nonlinearity an ultracold atomic wave packet may form a many-atom 

SIh , soliton guided by a travelling or standing laser mode of constant frequency. Thus the atoms are immune 

O ■ to wave packet spreading if no additional forces are present. Evidently the realization of such solitons is 

of great importance for atom interferometry and other domains of quantum optics. 

The considerations mentioned above are in some sense incomplete. Atoms are massive and therefore 
necessarily subject to the influence of the gravitational field of the earth. This is one of the main differences 
between nonlinear photon optics and nonlinear atom optics. In the set up of |6j it is the acceleration 
5^ \ due to gravitation which gives the outcoupled atoms a distinct direction. Will this destroy the soliton 

character of the matter wave packet, so that in fact the set up sketched above is of limited use? In any 
case a reshaping of the wavepacket is to be expected. Let us therefore look at the same problem from a 
different point of view: Consider an unaccelerated ultracold many-boson soliton (no force present) and 
study the influence of a frequency chirping of the laser wave. Will this allow us to manipulate the packet 
envelope so that a coherent control of the soliton shape becomes possible? It is the aim of this paper 
to show that this can indeed be done. To do so we have in a flrst step to reexamine all the underlying 
approximations used to derive soliton solutions in the force-free (or constant laser frequency) case. We 
will thereby assume the reader to be familiar with 

We discuss an ideal gas of identical two-level atoms interacting with the vacuum electromagnetic fleld 
and a travelling laser wave. The gas is assumed to be composed of ultracold Bose atoms, which can 
be treated as a vector quantum field with two components ijjg and ipe corresponding to the internal 
ground state and the excited state of the atoms, respectively. The internal energy difference is huja- Due 
to the coupling to the electromagnetic vacuum (virtual) photon exchanges between different ultracold 
atoms are possible. After the elemination of the vacuum field these photon exchanges result in a long 
range interatomic correlation which is responsible for the existence of atomic solitons. Following we 
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propose a scheme where ultracold atoms are loaded into a travelling laser beam with frequency w/ , which 
propagates in the z direction. But in contrast to the additional influence of a homogeneous force field 
parallel to the z direction is explicitly taken into account. This leads in the inertial reference frame to 
additional potential terms —mazipg/,, in the coupled differential equations for the ground and excited 
field components of the atomic gas whereby a may be positive or negative. 

Let us first note the equivalence of the influence of a homogeneous force field to a frequency chirping. 
The narrow atomic wave packet is assumed to propagate collinearly within the laser beam. The following 
transformations 

^gir,t) = 0,(r,i)e^'='^^e-*"(^'*\ (1) 

= n^+'>{x,y)e''"'-"^'* (3) 
where kc is the initial wave vector of the center of the atomic beam and where 

-M = -^^t-mazt+^^^'j (4) 

corresponds to the total kinetic energy of the atomic beam, yield the equations of motion for the atomic 
field envelopes (^g/e(r,f) and the Rabi frequency of the laser field 0(+^(r). These equations are of such a 
nature that after having changed to the co-moving frame of the accelerated atoms by 

z ^ z — fikct/m — ■^CLt^ ^ (5) 

the effect of the homogeneous force field is reduced to a chirping of the laser frequency. This means that 
after the transformations (|l|)-(|^) we obtain the same equations of motion for the envelope fields as in the 
forcefree case group velocity Vg = hkg/m = 0) but with a time dependent effective detuning A(t) 
due to the Doppler shifted laser detuning 

A{t) — uji — LOa ~ ujj. ~ kiiTikc I + at) (6) 

where ujr = Uk^ /2m is the recoil frequency. We introduce Aq = A(a = 0) for later use. 

To avoid the influence of spontaneous emission, which induces inelastic processes, we assume a large 
detuning. The probability of light scattering by single atoms, proportional to needs to be suf- 

ficiently small. However, for negative detuning the mean field interaction between atoms can be signif- 
icantly changed without excessive recoil or photoassociation losses due to the presence of vibrational 
quasimolecular states Q . 

In that region far from atomic resonance (that holds as long as |il/A(i)| ^ 1) we can adiabatically 
eliminate the excited state component up to the order A^^(t) and obtain the nonlinear Schrodinger 
equation of Ref. but with the time dependent detuning A(t). This adiabatic elemination is possible 
as long as the relative change of the detuning due to its time dependence is small, that applies for times 
t <^T with T = |Ao/fc;a|. The resulting equation of motion for the ground state field operator of the 
atoms contains an effective single particle potential V{x^y) ^ il{x,y)/A{t) (the usual dipole potential 
from the laser field) and a nonlinear part ~ A~^(t) due to the two-body dipole-dipole interaction. 

Essentially all N ultracold Bose atoms occupy the same state. This justifies the time-dependent Hartree 
approximation for the coherent evolution of the ultracold atoms The iV-atom wavefunction 

$(ri, r2, . . . , r^v; t) is written as a product of N single atom wavefunctions $(ri, t) 

N 

$(ri,r2,...,rjv;i)~n$(r„i). (7) 

1=1 

To be consistent with the Hartree approximation (0) we neglect dissipative processes due to sponta- 
neous emissions and laser driven collisions between the ultracold atoms. Spontaneous emission can be 
ignored for a highly detuned laser (|A(t)| 3> 7, where 7 is the spontaneous emission rate) as shown 
numerically in for a comparable situation, whereas the effect of binary collisions depends strongly 
on the laser detuning flS] . For a red detuned laser we must assure that the laser does not excite a bound 
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state level v, which can be built up due to the dipole-dipole potential between a ground and an excited 
atom. If the detuning of the laser is choosen between two vibrational levels with energy difference hSi, , 
this condition is fulfilled for t <Ti, with T^, = \5y/kia\ if the linewidths of the excited bound states 7^ are 
small compared with (5^ (^ivl^v <C 1). Otherwise the laser driven collision rate is increased, the colliding 
pairs get a higher momentum than the other atoms and cannot be described by the Hartree ansatz. 

Furthermore we replace the true dipole-dipole interaction by a contact potential pO[ | . This is not allowed 
in general and has been criticised ||l6[, but should be reasonable for ultracold atoms. Since the range of 
the dipole-dipole interaction is of the order of an optical wavelength A = 27r/fc;, the description of the 
full potential by its first moment is expected to be valid for the long thermal de Broglie wavelenth of the 
ultracold atoms \dB ^ In practice the spatial orientation of the atomic dipole moment d is uncertain, 
so that we must take the mean value over the orientation of the dipoles in the dipole-dipole interaction. 
With the averaged dipole-dipole interaction we get a positive correlation volume Vc = 2n/Sk^. 

Thereby we obtain as an intermediate result a nonlinear Schrodinger equation for the effective single 
particle states of the ground state atoms: 

' dt 2m ^ 4A(0 2A2(t) 

The possibility of manipulationg the value and the sign of the scattering length, corresponding to the 
mean field interatomic interaction, is also studied in [l2| , where the neglect of the dissipativc terms and 
the contact potential approximation are justified in another approach. 

Let us now try to find a solution of the nonlinear Schrodinger equation with a propagation of the atoms 
guided by the laser beam. Since the single particle potential in ^ is independent of the longitudinal 
coordinate z, we split the wave function in a transverse and a longitudinal part by a separation ansatz 
$(r, t) = u{x, y, t)ip{z, t). The transverse motion of the single particle wave function is then determined 

by 

du{x^ = -^.(x y t) - y t) (9) 

It contains a potential term which depends on the transverse laser profile fl{x,y) and decreases with 
increasing detuning. For an appropriate profile, which may for example be a Gaussian, one can realize 
for the transverse wave function bounded energy eigenmodes uo(x, y, t). We assume the atomic ensemble 
to be in one of these eigenmodes. Due to the adiabatic theorem it will stay in this eigenmode collecting 
only an additional dynamical phase factor which will not be relevant as long as eT ^ h is fuUfilled 
whereby e is the energy difference to the next energy eigenstate and T is the typical time scale for the 
time dependence of the single particle potential introduced above. Turning to the particular case of an 
external force acting, it is worth to note that if the force ma is not collinear to the laser momentum then 
respective components in x or y direction result only in a displacement of the atomic center of mass in 
these directions, which can be compensated by an appropriate shape of the laser profile. 

In addition to this transversal confinement it is the longitudinal behaviour which is of central impor- 
tance for the shaping of ultracold atomic solitons and their stability in external force fields. We then 
finally obtain the following nonlinear Schrodinger equation in the co-moving accelerated frame 

where 

<t) = JdxJ dy\u,{x,y,t)\^\n{x,y)\' (11) 

is the nonlinearity coefficient which has become time dependent because of the influence of the force 
field. It corresponds to the Kerr-type nonlinear susceptibility in nonlinear optics. The coefficient demon- 
strates how with increasing effective detuning |A(t)| the dipole-dipole interaction is switched off and the 
nonlinearity in ( |lo|) disappears. 

The important consequence of acceleration or chirping is that the nonlinearity coefficient k is now time 
dependent. An adiabatic treatment is a priori not justified in a nonlinear equation. Therefore it is very 



3 



useful to make the following transformation to separate the time dependence from the nonlinearity 



2m,. 



K{t)N\ 



-1/2 



(12) 



We obtain 



. 2m dip 



+ i r(i)(^ = 0. 

n 



9^ 

This reveals the underlying physics: The time dependence of k, gives rise to an additional damping 
coefficient 

n{t) 

which is of the order T^^ . An increasing nonlinearity coefficient K[t) reflects according to (^) a growing 
longitudinal confinement due to the laser induced nonlinearity, and leads to antidamping. This is the 
case for the red detuned laser field propagating in reverse direction of the force field (a < 0). If both 
directions agree (a > 0) damping occurs. 



(13) 



(14) 



For times t with |r|< <C 1 that is for t <t^T the damping term in (13) may be treated perturbationally 
prf . Returning to the longitudinal part of the wave packet envelope we get (for a ^ still with reference 
to the co-moving accelerated frame) a normalized soliton solutions of the following form: 



mK(t)N , ,mK(t)N , 

sechi — — z \ exp 



2h 



dt' 



(15) 



with Ko = K{t = 0). This is the intended result. Under the influence of a homogeneous force field 
acting on the atoms or for a linearly chirped laser, the soliton changes its shape adiabatically in time 
due to the time dependence of the nonlinearity coefficient K{t). This opens the possibility of a coherent 
control of the shape of the soliton. The relative change of the width W{t) = \/2h/mK{t)N is constant 
W/W = — 2r~^sign(aAo). The related relative change of the height of the soliton is given by the inverse 
expression. The damping and antidamping mentioned above can be read off directly. If in the case of 
acceleration we go back to the inertial reference frame we see that the pulse maximum is also accelerated 
due to the homogeneous force. 

Let us now turn to the experimental realization and estimate the influence of the earth gravitational 
field. We consider the [351/2 ~* ■^Pi/2] transition of sodium atoms with wavelength 589nm {ki — lO^m"^) 
and a natural linewidth of 7/27r — lOMHz. We assume a red detuning Aq of lOOGHz and for the Rabi 
frequency flo ~ Aq/IO so that the adiabatic elemination of the excited atom state is justified. To keep 
the relative change of the detuning small and to justify our approximations, we have to restrict to times 



t smaller than T — \A/kia 
states and Si, ~ lOGHz [n5 



10 s. In fact, with 7^ = 20MHz as the linewidth of molecular bound 



the cold collisions mentioned above can only be avoided during a smaller 
time up to approximately T^, — \6^/kia\ ~ 100s since "f^/Si, <C 1, which is still a rather long travelling 
time for the atomic sample. To discuss the transversal confinement, we assume a Gaussian laser profile 
il{x, y) = r^o exp[— (x^ + y'^)/2bf] with a width bi of the laser beam much larger than the width of the 
atomic beam. Accordingly we obtain approximately harmonic oscillator eigenfunctions for the transverse 
atomic mode uo{x,y,t) with a width ba of the fundamental eigenmode. It can then be shown that 
the condition eT ^ h necessary for the adiabatic treatment of the transversal modes amounts to the 
condition |fc;a/rig| <C {ba/bif' <^ 1 which can always be fulfilled. For the experimental specification given 
above all approximations are therefore justified and our result (|l^) may be applied. To demonstrate how 
effective the confining nonlinearity still is, let us assume that the influence of the gravitational force field 
maintains for = 100s. The width and the height of the soliton is then changed only by 20%. This 
is compared with the spreading of an ordinary Schrodinger wave packet with inital width Wq under 
the same conditions: W{t) = (1 + {uJrtfikiWay^Y/'^. For Woh ~ 100, which is consistent with the 
initial width for the soliton in our case {N ^ 10^), the width of the Gaussian increases by about 40% in 
4 * 10~^s. In the time Ti, = 100s the width has increased by a factor of 10*. 

We turn to the positive aspect of an active manipulation of the soliton shape. For this the effective 
detuning may be changed by chirping of the laser frequency. If the resulting time dependence is linear, 
it can be described by an effective acceleration a. Then, to obtain a change of the soliton width of 20% 
in already Is of chirping, an effective acceleration of a = lOOg is necessary. 



4 



We thank F. Burgbacher and K.-P. Marzlin for fruitful discussions. M.H. also acknowledges helpful 
discussions with Y. Castin and financial support by the Optik Zentrum Konstanz. 



References 

[1] M.H. Anderson et al, Science 269, 198 (1995). 

[2] K.B. Davis et al, Phys. Rev. Lett. 75, 3969 (1995). 

[3] M.O. Mewes et al, Phys. Rev. Lett. 77, 416 (1996). 

[4] C.C. Bradley, C.A. Sackett, and R.G. Hulet, Phys. Rev. Lett. 78, 985 (1997). 
[5] R.J. Ballagh, K. Burnett, and T.F. Scott, Phys. Rev. Lett. 78, 1607 (1997). 
[6] M.O. Mewes et al, Phys. Rev. Lett. 78, 582 (1997). 

[7] M. Lewcnstein, L.You, J. Cooper, and K. Burnett, Phys. Rev. A50, 2207 (1994). 

[8] W. Zhang, and D.F. Walls, Phys. Rev. A49, 3799 (1994). 

[9] W. Zhang, D.F. Walls, and B.C. Sanders, Phys. Rev. Lett. 71, 60 (1994). 
[10] G. Lenz, P. Meystre, and E.M. Wright, Phys. Rev. A50, 1681 (1994). 
[11] G. Lenz, P. Meystre, and E.M. Wright, Phys. Rev. Lett. 71, 3271 (1993). 

[12] P.O. Fedichev, Yu. Kagan, G.V. Shlyapnikov, and T.M. Walraven, Phys. Rev. Lett. 77, 2913 (1996). 
[13] E.M. Wright, Phys. Rev. A43, 1066 (1974). 

[14] K.J. Schernthanner, G. Lenz, and P. Meystre, Phys. Rev. A50, 4170 (1994). 
[15] K. Burnett, P.S. Julienne, K.-A. Suominen, Phys. Rev. Lett. 77, 1416 (1996). 
[16] O. Morice, Y. Castin, and J. Dalibard, Phys. Rev. A51, 3896 (1995). 

[17] Y. Kodama, and M.J. Ablowitz, Stud. Appl. Math. 64, 225 (1980); Y. Kodama, and A. Hasegawa, 
IEEE J. Quant. Elec. QE 23, 510 (1987). 



5 



